Modes of counterion density-fluctuations and counterion-mediated attractions 

between like-charged fluid membranes 
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Counterion-mediated attractions between like-charged fluid membranes are long-ranged and non- 
pairwise additive at high temperatures. At zero temperature, however, they are pairwise additive 
and decay exponentially with the membrane separation. Here we show that the nature of these 
attractions is determined by the dominant modes of fluctuations in the density of counterions. 
While the non-pairwise additive interactions arise from long-wavelength fluctuations and vanish 
at zero temperature, the short-ranged pairwise additive interactions arise from short-wavelength 
fluctuations and are stronger at low temperatures. 
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Counterion-mediated attractions play a significant role 
in many physical and biological phenomena JT|— [To[] . These 
attractions are responsible for fitting DNA inside a small 
biological container such as a viral capsid or a nuclear 
envelope Pill]] and can also be crucial in promoting ad- 
hesion and fusion of biological membranes [0 . Accord- 
ingly, significant effort has been expended in develop- 
ing a practical way of investigating the nature of these 
counterion-mediated attractions. Besides an integral 
equation method [ jT3| , there have emerged two distinct 
approaches. The first approach ]i4| , |l5j ], based on a charge 
fluctuation picture, suggests that these attractions are 
mediated by correlated fluctuations of ion clouds of coun- 
tcrions. This approach is consistent with our conception 
of counterions as fluctuating objects and thus merits sig- 
nificant consideration |l|,||JJJ^,P • In the second approach, 
based on a zero temperature picture |p^pO| , the appear- 
ance of attractions between like-charged molecules is at- 
tributed to the strong charge correlations that drive the 
systems, together with counterions, into an ionic crystal. 
At first glance these two approaches appear to be con- 
tradictory to one another, but there is an evidence that 
they can, in fact, be complimentary [^l|-p3|. Despite this, 
there still remain fundamental discrepancies between the 
two that have yet to be resolved. For the case of two pla- 
nar surfaces a distance h apart, the zero-temperature pic- 
ture leads to an attractive force that decays expenentially 
with h In the charge-fluctuation approach, however, 
the attractive force decays algebraically as h" 3 as 
long as h is sufficiently large. Most recently, the appar- 
ent discrepancy between these two approaches has been 
examined [pij by Lau et al. Lau et al. have shown that 
the zero-temperature quantum fluctuations give rise to a 
long-ranged attraction which varies as hr 1 1 2 . Their the- 
ory, however, does not explain the crossover between the 
long-ranged interactions driven by thermal fluctuations 
and the exponentially-decaying interactions expected at 
low temperature. At finite temperatures, the long- 
ranged interactions supported by the charge-fluctuation 
approach should exist and constitute a dominant con- 



tribution to the plate-plate interactions. Furthermore, 
when applied to many rod systems p[ p^p5[ |, the charge- 
fluctuation approach suggests that interactions between 
rods are not pairwise additive while the exponentially 
decaying interactions between plates as implied by the 
zero-temperature approaches are pairwise additive. 

In this Letter, we present a simple theoretical approach 
to counterion-mediated interactions between fluid mem- 
branes in order to bridge the gap between the two ex- 
isting theories. In particular, we show that the nature 
of these attractions is controlled by the dominant modes 
of fluctuations in the density of counterions. At high 
temperature, the membrane interactions are dominated 
by Zong-wavelength fluctuations in the counterion den- 
sity, i.e., fluctuations at large length scales. The result- 
ing Zong-wavelength interactions are shown to be long- 
ranged and non-pairwise additive. As the temperature 
decreases, the high-temperature behavior of the mem- 
brane interactions crosses over to the behavior expected 
at low temperatures. At low temperatures, the mem- 
brane interactions are mainly determined by the short- 
wavelength charge fluctuations. We find that the result- 
ing interactions decay exponentially with the separation 
between the membranes and are approximately pairwise 
additive. Finally, we obtain a phase diagram to depict 
the two distinctive regimes characterized by the corre- 
sponding dominant modes, and the crossover boundaries 
between the two. 

The system we consider here consists of negatively 
charged parallel membranes surrounded by neutralizing 
counterions of opposite charge. In the following deriva- 
tion, we assume that counterions are localized in the 
plane of the membrane. Thus our approach presented 
here is relevant to the strongly charged case. More pre- 
cisely, the separation between plates is much larger than 
the Gouy-Chapman length, a typical length scale within 
which counterions are confined. The main purpose of 
the present work is to study the crossover from the 
high-temperature results for the membrane attractions 
to the behavior expected at low temperatures. Since this 
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crossover occurs at low temperaures or at high densities 
of counterions, the assumption of localized counterions is 
reasonable. 

The charge distribution on a layer j is described by 
the local surface charge density, &j(r±) = —eao + emZ, 
where e is the electronic charge, <tq is the average coun- 
terion number density, m = 0, 1, 2, 3, etc., is the number 
of counterions per unit area at rj_ = (x, y), and Z is the 
counterion valency. In order to calculate the free energy 
and the corresponding charge correlations (cf. Eq. |^), we 
use two-dimensional Debye-Hiickel (DH) theory. This ap- 
proximation, however, fails to capture the strong charge- 
correlations at low temperatures. This defect in the DH 
theory has been corrected in an approximate way by tak- 
ing into account short-ranged charge correlations over 
the size of ions 1 21 - 23 1 . We thus implement DH theory 
with the counterion size via the two-dimensional form 
factor g(r_L, r' ± ) = 6(|rj_ — r^_| — D) /ttD 2 where D is the 
diameter of the counterions. We find that the charge- 
fluctuation contribution to the free energy per area is 
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where A -1 = 2ttZ£b<Jq is the inverse Gouy-Chapman 
length and Ib = e 2 /eksT is the Bjerrum length, i.e, the 
length scale at which the electrostatic energy between 
two charges is comparable to the thermal energy. The 
matrix Q(kjJ is defined by the matrix elements 



f(k±) = ^. 



ffu( k J 
Afc, 



, — k I hi 



(2) 



where <?y(k_i_) is g(kj_) = 2Jl fc ^' P ^ if i = j and 1 other- 
wise, and J\(x) is the first-order Bessel function of the 
first kind. Finally, hij is the separation between plates i 
and j 

First note that the free energy of an A-plate system is 
not simply a pairwise sum of the corresponding two- plate 
results over all pairs of plates. Thus the pairwise additiv- 
ity is not always satisfied as will be detailed later. For the 
N = 2 case, our result in Eq. [I] reduces to the previous 
two-plate result [see Eq. 3 of Ref. ||], if D is set to zero. 
When D = 0, the free energy has a single minimum at a 
nonzero value k± — k< <C lA -1 for all values of A. The 
dominance of the long-wavelength charge fluctuations is 
responsible for the breakdown of the pairwise additivity 
of electrostatic interactions between macroions. It has 
been shown that the pairwise additivity for the case of 
charged rods breaks down if the expansion of the cor- 
responding interaction free energy in powers of £b di- 
verges ^|,p||22]|. For rod systems, the free energy is dom- 
inated by the zero-A: mode and thus this expansion con- 
verges only when the charge fluctuation along the rods is 
vanishingly small. For the two-dimensional case, the con- 
vergence of the ^-expansion can be tested by estimating 
6 = (fc^) _1 A _1 ; the ^-expansion is convergent if S < 1. 
We, however, find that 8 is smaller for smaller values of 



A -1 and is comparable to unity if A -1 < 1CP 3 . As in the 
rodlike systems, the pairwise additivity is violated unless 
the inverse of the Gouy-Chapman length is very small. 
Thus the pairwise additivity can easily be violated even 
in two-dimensional systems. 
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FIG. 1. The wavevector dependence of the interaction free 
energy between two plates separated by h — 5 A, as a function 
of k±. We have chosen D = 5A. When A" 1 = lA^ 1 , the free 
energy has only one minimum at k± -C lA -1 , while when 
A" 1 = 5.9A -1 the free energy has another local minima at 
k ± = k> = 0(lA _1 ). When A -1 = 7A" 1 , the second minima 
at large k± is overwhelmingly dominant. When A -1 = 10A _1 
the free energy has two local mimina at k± = 0(lA - ). 

The above analysis, however, cannot be pursued when 
the temperature is low or when the density of ions is high. 
In this case, it is crucial to incorporate D ^ [|l] p3[ ; 
the ion size constitutes an important length scale at low 
temperatures or at high density of ions. Having under- 
stood that the dominance of long- wavelength charge fluc- 
tuations for the case D = leads to the breakdown of 
pairwise additivity, we now examine the low-T behavior 
of the free energy (thus with D set to a finite value) in 
Fourier space. Let us first consider the case of two plates 
seperated by a distance h. For convenience, we con- 
sider the following quantity: Q(kj_) = k± log[det Q(k±)], 
i.e., the first term in {. . .} of Eq. ^ mutiplied by k±. 
In Fig. [l], we plot this quantity Q(k±) as a function 
of k± for several different values of the Gouy-Chapman 
length A. We have chosen D = 5A and h = 5A. 
When A -1 = lA _1 , Q(kj_) has a single minimum at 
k± = fc< <C 1. This implies that the free energy is dom- 
inated by long-wavelength charge fluctuations as in the 
previous case of point charges. We find that the func- 
tion Q(kj_) has two minima at k± = k< <c lA^ 1 and 
at k± = k> = OilA- 1 ), respectively, for A" 1 > 4.2A" 1 
(not shown in the figure). As A -1 changes, the minimum 
at k^_ ^— k varies monotonically and is deeper for larger 
A -1 . The minimum at kj_ = k±, however, is roughly 
independent of A. When A -1 w 5.9A -1 , the two min- 
ima are comparable in magnitude. When A -1 = 7A _1 , 



2 



the function Q(kj_) is overwhelmingly dominated by the 
second minimum at k± = kj_, as shown in the figure. 
At A -1 = A^- 1 = 7 '.2A~ X , the second minimum diverges. 
This is suggestive of the onset of crystallization of coun- 
terions as will be detailed later. Even though the region 
A -1 > A^- 1 is certainly beyond the validity of our the- 
ory it is neverthless interesting to see what our theory 
impiles for that region. Notably, the free energy curve 
corresponding to A -1 = KM -1 has two local minima at 
large k± = The existence of multiple minima 

at large k± assures that the system is in a solidlike phase. 

Our results in Fig. [l] imply that there are two distinct 
contributions to the free energy: long-wavelength fluctu- 
ations and short- wavelength fluctuations in the density of 
counterions. They also imply that the short-wavelength 
fluctuation contribution to the free energy has a much 
narrower peak if A -1 3> 5.9A -1 . This enables us to sep- 
arate the short-wavelength contribution from the long- 
wavelength contribution. By noting that e~ k±h does not 
change appreciably over the region inside the peak at 
k± = we find, up to /i-independent terms 
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where C,(x) is the zeta function (thus £(3)/167r ~ 0.024). 
The first term denoted by Fjsw is the free energy cal- 
culated with D set to zero |?]-||,|27j and is the long- 
wavelength free energy. Our previous analysis on Q(kj_) 
implies that the short-wavelength free energy denoted by 
Fsw, i.e., the second term in Eq. ||, is dominant over the 
long-wavelength free energy Flw at low temperatures 
(A -1 S> 5.9A -1 ) and decays exponentially in space. The 
exponentially decaying interaction between two pla tes is 
consistent with the zero-T analysis in Refs. [l^]24] |. At 
high temperatures corresponding to A -1 <C 5. 9A^\ how- 
ever, the free energy is mainly determined by Flw In 
this case, the fluctuation contribution to the pressure be- 
tween the two plates scales as — /i~ 3 ]7|-p|,p7[| . Also note 
that this long- wavelength contribution vanishes at T = 0. 
This follows from the fact that Q(k±) is roughly inde- 
pendent of A. To estimate the temperature dependence 
of Fsw, note that the prefactor of this term varies as 
T -1 . The fc_L-integral of this term depends on the depth 
and width of the second minimum of Q(k±) at k± = k^_. 
While the width is roughly independent of A for given h, 
the minimum becomes deeper with the increasing A -1 (or 
decreasing T) . This proves that Fsw is more negative at 
low T than at high T, as opposed to Flw- 

For a given value of h, there exists a special value X cr 
at which the crossover between the two distinctive be- 
haviors of the plate interaction (Flw and Fsw) takes 
place. By requiring Jj^ (Flw — Fsw) = 0, we have the 
following transcendental equation for \ cr : 
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To solve the trascendental equation, we have chosen 
D = dA. Fig. U describes distinct regimes characterized 
by the corresponding dominant modes of fluctuations, 
and the crossover boundaries between them; the regimes 
where the long-wavelength (LW) and short-wavelength 
(SW) fluctuations dominate are denoted by LW and SW, 
respectively. When A -1 is smaller than 2.7A~ 1 , the plate 
interaction is solely determined by the LW fluctuations 
for the whole range of h. At A" 1 ~ 2.7A~ 1 , marked 
by the vertical dotted line on the left, the SW fluctua- 
tions start to contribute to the plate interaction. When 
2.7A~ 1 < A -1 < 7.6A -1 , however, the plate interaction 
is determined by the competition between the two; the 
crossover from the LW to SW regime takes places for 
larger value of h at low temperatures (corresponding to 



larger A -1 ). At A" 



a; 



7.6A" 1 EM, marked by 



the vertical dotted line on the right, the crossover occurs 
only when h — > oo. Beyond A^ 1 , the SW fluctuations 
solely determine the plate interaction. 
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FIG. 2. Phase diagram for two plates. The regimes where 
long-wavelength and short-wavelength fluctuations dominate 
are denoted by LW and SW, respectively. 

The appearance of two distinctive competing interac- 
tions, i.e., Flw and Fsw, can also be understood in 
terms of in-plane charge correlations for a single plate: 
Gi(r X) r^) = (a(r ± )a(r' ± ))-{a(r ± )) (a(r' x )). The long- 
wavelength contribution to the charge correlation was 
shown to scale as Glw( v X, r '±) ~ — ksT/\rx — r^J 3 f p7| . 
Note that this correlation vanishes as T — > 0. Our re- 
sult in Eq. |l| implies that the short-wavelength charge 
correlation function is given by 
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where Jo(x) is the zeroth-order Bessel function of the first 
kind. Unlike the long-wavelength correlation, the behav- 
ior of the short-wavelength correlation is determined by 
the nature of the poles of [f + Xk±/g(k±)] . We find 
that Gsw shows an oscillatory decay, and that the am- 
plitude of Gsw varies as T- 1 [1 + (fc>A)- 1 . 9 (fc>)] and 
is larger at low temperatures. For the two-plate case, this 
oscillatory charge correlation essentially gives rise to the 
SW interaction. 

At high temperatures, we can consider each plate 
to consist of large domains, i.e., counterion-rich and 
counterion-poor domains. The size of the domains is 
on the order of 3> lA, and thus these domains 

can form huge dipoles, resulting in a long-ranged attrac- 
tion. The long- wavelength fluctuations couple over many 
plates, leading to breakdown of pairwise additivity pSy 
in the TV-plate case. We also find that the power-law be- 
havior of the in-plane charge correlation crosses over to 
an exponential decaying form as N — » oo |2{|. This is 
clearly suggestive of the breakdown of pairwise additivity. 
As the temperature decreases, however, the plate interac- 
tion driven by long-wavelength fluctuations crosses over 
to a distinct behaviour controlled by non-monotonically 
decaying charge correlations. At low temperatures, each 
domain becomes overall charge neutral and thus the dis- 
tiction between domains is meaningless. In this case, the 
local correlation between a counterion and a backbone 
charge in its neighborhood dominates the free energy. 
There is thus strong cancellation of repulsions (between 
like charges) with attractions (between opposite charges). 
This results in an exponentially decaying, short-ranged 
attraction between the plates. 

In conclusion, we have presented a systematic ap- 
proach to study the nature of counterion-mediated at- 
tractions between fluid membranes. We have shown that 
the nature of these attractions is determined by the dom- 
inant modes of fluctuations in the density of countcri- 
ons. At high temperatures, fluctuations at large length 
scales determine the membrane interactions; the result- 
ing interactions are long-ranged and not pairwise addi- 
tive. Charge densities of biomembranes range from 0.03 
to 0.24/nm 2 , corresponding to LW regimes at room tem- 
perature. It is thus clear that many-body, non-pairwisc 
additive interactions operate in biomembrane systems at 
room temperature. At low temperatures, however, the 
membrane interactions are dominated by SW charge fluc- 
tuations and are exponentially decaying with the mem- 
brane spacing. In this case, only the nearest pairs of 
plates couple strongly with each other. Surprisingly, this 
implies that the pairwise additivity is restored at very 
low temperatures (the non-pairwise additive interaction 
becomes smaller, and eventually vanishes, as T —> 0). 
The approach presented here allows one to systematically 
study the crossover, as the temperature decreases, from 
the high-temperature, long-ranged attractions to the be- 
haviors expected at low temperatures. 
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